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SUMMARY 


This report presents a robust control design using strictly positive realness for second- 
order dynamic systems. The robust strictly positive real controller allows the system to be 
stabilized with only acceleration measurements. An important property of this design is that 
stabilization of the system is independent of the system parameters. The control design 
connects a virtual system to the given plant. The combined system is positive real regardless of 
system parameter uncertainty. Then any strictly positive real controllers can be used to achieve 
robust stability. A spring-mass system example and its computer simulations are presented to 
demonstrate this controller design. 

1. INTRODUCTION 

Positive real (PR) systems have many applications for shape and vibration control of 
large flexible structures. In most of those PR designs, the output of the plant is usually 
assumed to include velocity, and the sensors are assume to be collocated with the actuators. In 
[1], position and velocity feedback are used together to control large space structures, and the 
controllers are strictly positive real. PR feedback with velocity measurement is examined in [2] 
for the control of a flutter mode. [3] presents a robust multivariable control of structures using 
a passive controller in which the velocity sensors are collocated with the control actuators. 
Several passive control designs using acceleration, velocity and position measurements are 
presented in [4]. [5] generalizes the designs in [4] to handle nonlinear systems. The method 
presented in [6] uses displacement sensors. Similarly, [7] examines direct position plus 
velocity feedback. A feedforward positive real design can be seen from [11]. 

Nevertheless, in some application areas, only acceleration is directly measurable. Even 
though velocity and position may be obtained by integrating the measured acceleration, bias in 
velocity and position will decrease the accuracy of the integration. Therefore, in this study we 
develop a robust controller for multivariable second-order system when only acceleration is 
directly measurable. 

In this report, we review some definitions and a theorem associated with dissipativeness 
and passivity. Dissipativeness and passivity are then related to strictly positive realness and 
positive realness. Using these backgrounds we develop a virtual system to compute an output 
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which will make the combined system of the plant and the virtual system positive real (PR). 
The inputs to the virtual system are only acceleration and the control force applied to the plant. 
More important, the virtual system is model independent and thus the global system is robustly 
positive real. Therefore the input / output controller can be constructed by any strictly positive 
real controllers. When the stiffness matrix of the second-order system is positive definite, we 
show that it is possible to stabilize the displacement if the actuators are properly located. With 
this design, the displacement is globally asymptotically stable. A spring-mass example with 
three masses and no damping is used to illustrate our design method. Computer simulations are 
also presented. 

2. PRELIMINARIES 

The concept of dissipativeness describes an important input-ouput property of dynamical 
systems. Consider a system with input u and output y, where u is an mxl vector and y is a 
pxl vector. A supply rate for the system is defined as follows. 

Definition 1 [8]: A supply rate is a real function of u and y defined as 

w(u, y) = y T Q y +2 y T S u + u T R u (1) 

where Q, S, and R are constant real matrices with dimensions pxp, pxm and mxm 
respectively. 

Q and R are usually symmetric matrices. w(u,y) is often called the input energy into the 
system. Dissipativeness is defined with respect to the supply rate w(u,y) in the following 
definition. 

Definition 2 [8]: The system with input u and output y is called dissipative with respect to 
the supply rate w(u,y) if for all locally integrable u(t) and all T> t 0 , we have 

J w(t)dt£0 (2) 

l o 

where x(t 0 )=0 and w(t)=w(u(t),y(t)) is evaluated along the trajectory of the system 
interested. 
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Eq.(2) means that an initially unexcited system can only absorb energy as long as the 
system is dissipative. If the supply rate represents the input energy into the system, then Eq.(2) 
states that a system with no initially stored energy transforms the input energy into either stored 
energy or dissipated energy. Thus no energy can be generated from a dissipative system. 

Passivity is defined as a special case of dissipativeness. 

Definition 3 [8]: A system is passive if and only if it is dissipative with respect to the supply 
rate 


w(u, y) = u T y 


(3) 


An algebraic condition for passivity can be found if the system is represented by the 
state-space equations 


x = f(x)+G(x)u 

y = h(x) + J(x)u (4) 

where f(x) and h(x) are real vector functions of the state vector x, with f(0)=0, h(0)=0, and 
G(x) and J(x) are real matrix functions of x. These four functions are assumed to be infinitely 
differentiable. We also assume that u and y have the same dimension. The system is 
furthermore assumed to be completely controllable. Theorem 1 provides a test for the passivity 
of a system written in the form of Eq. (4). 

Theorem 1 [9]: The system is passive if and only if there exist real functions <J>(x), l(x) and 
W(x)with <(» (x ) continuous and with 

<|>(x) £0, for all x (5) 


and 


< K 0)=0 


(6) 



such that 


(i) V T <|>(x)f(x) = -l T (x)l(x) 

(ii) 1 / 2G T (x)V<{>(x)=h(x) - W T (x)l(x) ( 7 ) 

(iii) J(x) + J T (x) = W(x) T W(x) 

Moreover, if J (x) is a constant matrix, then W(x) may be taken to be constant. 

The function <j>(x) is generally not unique for a given dynamic passive system. 
Nevertheless, a physical meaning can be attached to it. It is shown in [9] that 

2 J u T (t)y(t)dt = <t>[x(T)] - <>[x(t 0 )] + J* [l(x) + W(x)uf [l(x) + W(x)u] dt (8) 

1 A ^ _ 


Eq.(8) may be interpreted as the conservation of energy equation, ^(x) is a stored energy for 
the system. The first integral corresponds to the input energy to the dynamic system. The 
second one is proportional to dissipated energy, and it is always nonnegative. As a 
consequence, Eq.(8) means that the energy input is equal to the variation of stored energy plus 
the loss of energy which is a positive function. 

A linear system is passive if and only if its transfer matrix is positive real [10]. Passivity 
can thus be seen as a generalization of positive realness for nonlinear systems. Since the 
systems investigated here are linear, we will equivalently use these two concepts for the rest of 
this report. 

3. A VIRTUAL SYSTEM DESIGN 

The multivariable system (Plant (P)) is described by 

Mx + Dx+Kx = Bu (9) 

where u is an mxl control vector, x is an nxl state vector, M is an nxn symmetric positive 
definite matrix, D and K are nxn symmetric positive semi-definite matrices, and B is an nxm 
matrix. Let a virtual system (V) be defined with the following equation 
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x =Ax +B’u 


( 10 ) 


where A is an Ixn matrix, B' is an lxm matrix, and x, is an lxl vector. The following 

Theorem 2 allows us to compute an output y that makes the global system (which is a 
combined system of the given plant and the virtual system) positive real. 


Theorem 2: Let H v , A and B be chosen such that 

2H v A = B t 
B t M^ =2 H v 

where M , is an 1 x 1 positive semi-definite matrix. If 

y = H v x. 


(ii) 


( 12 ) 


then the system with input u and output y is positive real. 
This scheme is illustrated in Fig. 1 


u 

Plant (P) 

Acceleration 

- . - 

System (V) 

y 
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Figure 1. A Virtual System 


Proof: For this proof, it is useful to represent the system with a state-space representation. Let 
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X T T T T T /T *T T* T \ 

= (Xj X 2 x' X 4 )= (x 1 X xi X a ) 

The equations describing the global system may be rewritten as 


( 13 ) 


jx = f(X) + G(X)u 
[Y = h(X) + J(X)u 


where 


f(X) = 


- m 'DXj- M -1 KXj 


A M ~ 1 D x 2 — AM "’Kx^ 


G(X) = 


0 

m'b 

0 

AM _1 B +B' 


h(X) = 


( o ^ 
0 
0 

H v x 4 y 


J(X)=0 


Let a candidate for the function <J>(X) in Theorem 2 be 


<t*(X) = yx T M x + ~x t K x + y (x, — Ax) T M, (x. - Ax) 


(15) 


(16) 


(17) 

(18) 

(19) 


where M , is positive semi-definite. The sum of the first two terms corresponds to the stored 
energy of the plant. The additional term is added for the positive real design. The function <XX) 
can be written using the state variables as 
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0(X)- 2 x 2 M x 2 + 2 x I Kx i + 2 ( x *~ Ax 2 ) T M.(x 4 - Ax 2 ) 


( 20 ) 


<KX) is a positive function and <{>(0)=0. It must be checked that there exists a function 1 (X) 
such that 


V T <KX)f(X)=-l T (X)l(X) 


This calculation is considerably simplified when we notice that 

vVx)f(X) = ^ 

u * 0 


( 21 ) 


( 22 ) 


As a consequence we have 


V>(X) f (X) = x T (M x+ K x) + j (x . - A x) T M. (x. - Ax) 

H-^x, -Ax) t M, (x.-Ax) 


(23) 


u=0 


When u = 0, the last two terms cancel out and therefore 
V T <(>(X) f(X) = x T (M x + Kx) I 

n = 0 

Thus we finally have 

V T <j>(X)f(X) =-x T Dx =-XjDx 2 


(24) 


(25) 


Since D is positive semi-definite, it is possible to find a matrix R such that D = R T R . The 
above equality becomes 


V T <t>(X) f(X) = - (Rx / (R x 2 ) = - 1(X) T 1(X) (26) 


where 1 (X)=Rx r Thus equality (i) from Theorem 1 is satisfied. Equality (iii) of Eq. (7) 
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reduces to 


J(X) + J T (X) = W T (X)W(X) = 0 (27) 


The function W(X) is therefore equal to zero. Equality (ii) of Eq. (7) becomes 

h(X) = ^G T (X)V<(>(X) 


(28) 


Only the partial derivatives with respect to velocity will be used to evaluate Eq. (28). We 

have 


d<j> T t 

^- = XjM + (AXj-x 4 ) M,A 

d$ T ( 29 > 

— = (x 4 -Ax 2 ) M. 


The function h(X) is such that 


2h(X) = (M 


B) 


d<}> 

8x 


a t 


+ (A M 


'b + b') t 


2 ) 


9 < 1 > 

dx 


\ T 


aJ 


(30) 


Obvious simplifications yield 

2h(X) = (B t -B T M^A)x 2 + B t M^x 4 (31) 


h(X) equals H v x 4 if the following equations are satisfied 

b t -b' t m^a = o 

B T mI=2H v 


( 32 ) 


Those equations can be rewritten as 
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( 33 ) 


2H v A = B t 
B T mI=2H v 


and the theorem is proved • 

There are several possible ways to solve the above system of equations. Given H v and 
B, we can solve for some possible A, M. and B\ At the end of the calculation, it must be 
checked that M, is positive semi-definite. Another method consists of choosing B, A and a 
positive semi-definite M , and then solving for possible B' and H v . 

4. CHOICE OF A CONTROLLER 

If the output of the global system is chosen as in Theorem 2, then the system is positive 
real. Thus the closed-loop system is uniformly asymptotically stable with zero input if the 
controller is strictly positive real [3]. That is, for this case, we have 

lim (H v x,) = 0 (34) 


Our goal is to let x go to zero. Theorem 3 may be used to achieve this goal. 

Theorem 3: Assume that Theorem 2 is used to make the global system PR. Furthermore 
assume that 

T 

(t) B x = 0 and u = 0 imply x = 0 . 

(ii) K is positive definite. 

(iii) The system is connected to a dissipative closed-loop controller. 

Then lim x(t) = 0. 

I — * *■ 


Fig.2 shows the control scheme for the plant (P). 
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Figure 2. A SPR Controller for the Virtual System 


Theorem 3 allows us to design a robust controller for Plant (P). No knowledge of the 
constant matrices M, D and K is required. Furthermore, the only measurements needed are 
acceleration and input. Acceleration may easily be measured for many practical systems by 
using common accelerometers. The input u may be obtained by measuring the output of the 
SPR controller. 

The proof of Theorem 3 uses the following two lemmas. 

Lemma 1: Assume that the Laplace transform of f(t) and f(t) exist in a neighborhood of the 
origin. Furthermore assume that lim f(t)=0.Then lim f(t) = 0. 

I — » •• t — * •* 

Proof: Let F(s) be the Laplace transform of f(t). The final value theorem yields 

lim sF(s) = 0 (35) 

»-» o 

The Laplace transform of the derivative of f(t) is 

L[f(t)] = sF(s)-f(0) (36) 

As a consequence we have 
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( 37 ) 


lim f(t) = lim s (sF(s) - f(0)) = 0 

t -+ ~ t -* 0 

Lemma 2: Let e(t) be a function of time and let e(t) go to zero as time increases. Then if 
satisfies the differential equation 

Dx + Kx =e (38) 

where D is positive semi-definite and K is positive definite, then x converges to zero. 

Proof: Let m denote the rank of D. There exists an invertible nxn matrix P such that 

D* = P DP 1 (39) 


where 


D* = 


0 0 ] 

o dJ 


(40) 


D 22 is an mxm positive definite matrix. Let K* be defined as 


K =P K P' 


(41) 


K* may be written as 


K = 


k„ k; 2 
k;, k* 


(42) 


The dynamical equation can be written as 


P D P" 1 (P x) + P K P" 1 (P x) = P £(t) 


( 43 ) 


Let y =Px and tl (t ) = P £(t ) . The system is now described by 



these equations are reduced to 


if y = 


V 

w 


and T) = 


'V 

L*n J 


K *u Yi + k* 12 y 2 = “n , (0 

1^22 +K 2»yi ^ 12 ^ 1 ~ ^2^^ 


( 45 ) 


The first equation can be solved in terms of yj and Eq. (45) reduces to 

, yi c - K M ,K i 2 y,+K^ 'n,(t) 

D„ y, + (K’j - k;,k;;' k; 2 ) y 2 = 11,(0 + k;X;‘ n , <•> 


• • • “ 1 • 

D 22 and (K 22 - K 2J K n K 12 ) are positive definite matrices. Thus y 2 may be considered as 

the output of a strictly stable system. The output of the strictly stable system converges to zero. 
The parameter y 2 will therefore go to zero. The first equality in Eq.(46) shows that y also 

goes to zero. Consequently, y converges to zero and so does x. 

Proof of Theorem 3: (-u ) is the output of the dissipative controller. A dissipative controller 
is always strictly stable. Knowing that y goes to zero, we can therefore conclude that u also 
goes to zero. Furthermore, we have 

2H v x 4 =2H v Ax+2H v B'u (47) 

by multiplying Eq. (10) with 2H V . Since 2H v A = B T , this equation may be rewritten as 

B T x =2H,x i - 2H,B'u (48) 

2 H v B 'u goes to zero as u goes to zero. Furthermore, we know that y= H v x t converges to 
zero as time increases. Using Lemma 1 allows us to say that H v x t also goes to zero if we 

T • 

assume that the Laplace transforms of x t and its derivative exist. As a consequence, B x 
goes to zero. The equations describing the system are linear and consequently continuous. 
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T 

Thus, if B x and u go to zero, x goes to zero according to assumption (i) in Theorem 3. The 
dynamics of the closed-loop system is now 

Dx+Kx=Bu-Mx =e(t) (49) 

where e(t) vanishes as time increases. Using Lemma 2 we conclude that x ( t ) goes to zero. 

5. EXAMPLES 

We study the simple example of a system with three masses, three springs and no 
dashpots. The example is shown in Fig. 3. 


xl 


x3 


/I kl k2 *3 


m3 


Figure 3. A Spring-Mass System 


This system needs to be stabilized as it is not naturally asymptotically stable. With no control 
and non-zero initial conditions, the three masses oscillate since there is no damping. The 
equations describing the system in Fig.3 are 


mjXi+(k, + k 2 )x, -k 2 x 2 = u, 


" m 2 X2 — k 2 x, + (k 2 + k 3 )x 2 - k 3 x 3 = u 2 
m 3 X3 — k 3 x 2 + k 3 x 3 = u 3 


( 50 ) 


The matrices M, D and K are 
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( 51 ) 


M = 


m, 

0 

0 


0 0 
m 2 0 
0 m 3 


D = 0 


( 52 ) 


k, + k 2 



-k 2 0 

k 2 + k 3 — k 3 

“ k 3 k 3 


( 53 ) 


M and K arc positive definite as long as none of the masses and the spring constants is equal to 
zero. Several possible controller designs can be used here. 


m=n=l=3 


There arc three control parameters here. A reasonable choice is 


B = 


1 0 0 
0 1 0 
0 0 1 


( 54 ) 


and the control vector u is defined by u T = (u, u 2 u 3 ). Obvious solutions to Eq.(ll) are 
given by 


A =1 


3x 3 


h. = 1b t 

B'=XB 



3x 3 


( 55 ) 


where X is an arbitrary strictly positive real number. As a consequence, the vector x, is 
generated by the differential equation 
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x. = x + XB u 


(56) 


All the assumptions of Theorem 3 are satisfied. The vector x may therefore be controlled with 
the help of any dissipative feedback controller. A simple choice consists of taking a constant 
controller, that is a controller with a transfer matrix of the form k I, where I is the identity 
matrix. 

The simulation is made on MATLAB. The following values are used in the simulation: 


m,=m 2 = m 3 =l (57) 

kj = 1 k 2 = 2 k 3 = 3 (58) 

The initial conditions were arbitrarily chosen to be 

x, = 5 x 2 = - 2 x 3 =9 (59) 

0 0 0 

x j =3 x 2 =5 x 3 = - 4 (60) 

0 * 0 0 


For the vector x , , we choose the simple initial conditions 

x, = 0 x, = 0 (61) 

The constant X is equal to 0.5. The gain of the feedback controller is k = 1 . The plot of the 
displacements is shown in Fig. 4. In the following plots, xj is indicated by — , X 2 is indicated 
by ••• and X 3 is indicated by — . 
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Figure 4 


The control goal is achieved. The three displacements vanish with time. 
Nevertheless, this design requires that a force be applied on each of the masses. It is possible 
to reduce the number of actuators with the following control design. 


m=l=2 


Here only two forces are be applied to the system. Thus there are three possible choices, 
depending on what masses the forces are applied. Let us start with 


B = 


1 

0 

0 


0 

1 

0 


(62) 


This choice means that the forces are applied on the masses m t and m 2 . The control vector u 
is such that u T = (u , u 2 ) . The vector x . is now a vector with dimension 2. Eq.(l 1) has the 
following obvious solution 
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A = B 

h -It 
n v 2 2x2 

B'=U 


2*2 


m. = £i !x2 


(63) 


where X is an arbitrary strictly positive real number and I denotes the identity matrix. Thus x , 
can be computed from the following differential equation. 

x t = x + X.u (64) 


and the output of the system is y = — x , . 

A dissipative controller must be chosen to control the system. Here again, a constant 
controller is a simple possible choice. Its transfer function is k I, where k is a positive constant. 

- ■■ T * * 

It remains to ensure that Bx = 0 and u = 0 imply x = 0. If B x =0 and u=0, then the 
dynamical equations of the system become 

(kj + k 2 )Xj + k 2 x 2 =0 

<-k 2 Xj+(k 2 + k 3 )x 2 -k 3 x 3 =0 ( 65 ) 

m 3 X 3 - k 3 x 2 + k 3 x 3 = 0 


By differentiating the second equation and solving for x , we have 


_ k 3 i (k 2 + k 3 ) 


X3 = i Xl + 

k 2 




X2 


( 66 ) 


As x ] and x 2 are both zero, x 3 is also equal to zero. Thus the above equations are reduced 

to Kx = 0. Since K is positive definite, this yields x = 0. All the assumptions of Theorem 3 are 
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satisfied and we are now assured that x will go to zero. 


The closed-loop system is simulated with the same parameter choice as before. The plot 
of the displacements can be seen in Fig. 5. Here again the stabilization goal is achieved since 
the three displacements vanish as time increases. 



Figure 5 


It is possible to stabilize this system with a different distribution of forces. For instance, 
let 


B = 


0 

1 

0 


0 

0 

1 


(67) 


Two forces are applied respectively on mass 2 and mass 3. With the same design as above, a 
controller for the system can be designed. The plot of the displacements is presented in Fig. 6 
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with the same initial conditions. 



Figure 6 


Finally a third possible choice is 



1 

0 

0 


0 

0 

1 


( 68 ) 


In this case some forces are applied to mass 1 and mass 3. The plot of the simulation can be 
seen in Fig. 7. 
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Figure 7 


In all these case, the system is stabilized with the help of only two actuators. 


m=l=l 


Here we design a control system with only one actuator. This actuator may be located on 
any of the three masses. Let us first apply a force on mass 1, i.e. the matrix B is 


B = 


1 

0 

0 


(69) 


Eq.(l 1) in Theorem 2 has the following obvious solution 
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( 70 ) 


A = B T 

H V = 2 

B’= X 


M = 


1 

X 


where X is an arbitrary strictly positive real number. The state x, is calculated by integrating 
the differential equation 


x .= x , + 2 U 


(71) 


The output of the system is y = -j x , . 

Here again the SPR controller is chosen to be constant. Its transfer matrix is of the form 
G(s) = k, where k is any strictly positive real number. With this choice we are assured that x 
converges to zero. 


x . - 

It should be checked as before that B x=0and u = 0 imply x = 0 . The procedure is 
unchanged and once again those assumptions yield K x = 0. Since K is assumed to be positive 
semi-definite, x is necessarily equal to zero. 

The simulation is run with the same choice of initial conditions. The constant X is still 
equal to 0.5, and k is equal to 1. The three displacements go to zero as expected which can be 
seen in Fig. 8. 
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Figure 8 


The force could be applied on mass 3. The matrix B for this situation is 


B = 


0 

0 

1 


(72) 


The design method is unchanged. The closed-loop system has been simulated in this case with 
the same initial conditions and the same choice for the parameters involved. The plot is shown 
in Fig. 9. 
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Figure 9 


We could choose to apply the force on mass 2. In this case, 

B = 

Nevertheless, it can be checked that condition (i) of Theorem 3 is not satisfied in this case. 
Thus no controller design can be implemented with the above choice. 

6. CONCLUSIONS 

The control method presented in this report is particularly of interest for practical reasons. 
Only acceleration at certain locations of the system need to be measured by using common 
accelerometers. Furthermore, the design is model independent and no knowledge of the 
constants of the dynamic system is required. Finally, any strictly positive real controller can be 
used. Thus it is possible to choose one that yields a satisfactory transient response. 


0 

1 

0 
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